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Abstract. By accumulating all available radial velocity 
data and combining it with our own spectroscopy we pres- 
ent a new and improved ephemeris for Cygnus X-l (V1357 
Cygni/HDE226868). We initially calculate an orbital pe- 
riod of 5.599847±0.000018days and then, by including the 
more recent radial velocities of LaSala et al. (1998) and 
Sowers et al.(1998) we refine this to 5.599829 ± 0.000016 
days. Using 27 years' worth of photometry, we also calcu- 
late a photometric period of 5.599836 ± 0.000024 — re- 
markably similar to the spectroscopic one. 
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1. Introduction 

Cygnus X-l has become perhaps the most extensively 
studied X-ray binary, particularly due to its providing 
\ the first observational evidence for black holes. Webster 
■ & Murdin (1972) and Bolton (1972) independently associ- 
ated the X-ray source with HDE 226868, found an orbital 
period of 5.6 days and calculated the first mass function 
to indicate a compact object more massive than 3 Mq, 
the generally accepted maximum mass of a neutron star. 
The values of the masses of the two components have been 
refined a number of times, most recently by Herrero et al. 
(1995) giving 17.8 Mq for the supergiant companion star 
and 10.1 Mq for the black hole. 

A number of periods have been reported in the various 
lightcurves of Cygnus X-l. These include 294 days in the 
X-rays (Priedhorsky, Terrell & Holt 1983) and optical pho- 
tometry (Kemp et al. 1987), 150 days in the radio (Pooley, 
Fender & Brocksopp 1998) 78/39 days in UV polarimetry 
(Kemp, Herman, Rudy & Barbour 1977), optical photom- 
etry (Kemp, Herman & Barbour 1978) and 4.5 years in 
optical photometry (Wilson & Fox 1981). 

However, the orbital period is the only one that is 
confirmed in all wavebands and by most authors. It has 
been detected in the X-ray data of Ariel V/ASM (Holt, 
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Kaluzienski, Boldt & Serlemitsos 1976), RXTE/ ASM 
(Zhang, Robinson & Cui 1996) and CGi?0/BATSE (Pa- 
ciesas et al. 1997) and also in the radio (Pooley, Fender 
& Brocksopp 1998), optical photometry (Volushina, Lyuty 
& Tarasov 1997), optical spectroscopy (e.g. Gies & Bolton 
1982) and infrared (Leahy & Ananth 1992; Nadzhip et al. 
1996). 

The spectroscopic ephemeris has been revised many 
times using a variety of lines. The results of Gies & Bolton 
(1982) have been used (P orb = 5.59974 ± 0.0008 days) 
as the definitive values for over a decade and there has 
been little disagreement; while a variable orbital period 
was tentatively suggested (Ninkov, Walker & Yang 1987) 
this has never been confirmed and appears unlikely. It is 
only recently that a new ephemeris has been calculated 
(Porb = 5.5998 ±0.0001 days) and this is within the errors 
of the Gies & Bolton (1982) result (LaSala et al. 1998). A 
more accurate result is given as P or b — 5.59977 ± 0.00002 
days (Sowers et al. 1998) which also takes Hipparcos pho- 
tometrical data into account. 

The spectrum is that of a typical 09 supergiant 
with moderate Ha emission, presumably coming from 
the supergiant (e.g. Brucato & Zappala 1974), and also 
He II A4686 emission. Both of these emission lines are su- 
perimposed on absorption components and, by distin- 
guishing the two components, it has been shown by a 
number of authors (e.g. Aab 1983a) that the Hen emis- 
sion component is ~ 1 15 — 120° out of phase with the He II 
absorption component and the supergiant. As a result of 
this it is thought that the He n emission originates in the 
accretion stream. 

Neither of these absorption components are suitable 
for measuring radial velocity due to the inaccuracies in- 
volved in the removal of the superimposed emission. In- 
stead, as previously, we turn to the Hei absorption lines 
which are uncontaminated by emission. 

2. Observations and Discussion 

2.1. Spectroscopy 

Our spectra were obtained in June/ July 1997 using the 
Coude spectrograph of the Crimean Astrophysical Obser- 
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vatory's 2.6 metre telescope. The detector was a CDS9000 
(1024x256 pixels) CCD array. All observations were made 
in the second order of a diffraction grating with reciprocal 
dispersion of 3A/mm and resolution of 25000. The typ- 
ical exposure time for each spectrum totalled 1.5 hours 
resulting in a S/N of ~100. 

With a spectral width of 60A, our 20 spectra centred 
on He II A4686 also included the He I A4713 line. This is the 
line we have used for our radial velocity studies, although 
other sources have used alternative helium lines. 

The spectra were reduced using standard flat -field nor- 
malisation and sky subtraction techniques. Wavelength 
calibration was achieved using ThAr comparison spectra 
and to an accuracy of less than 0.5 km/s. The radial ve- 
locities were calculated by fitting a Gaussian to the core 
of each line and values of V r can be found in Table 1. 

To obtain the orbital solution we initially considered 
all of the radial velocity measurements. It was necessary 
to subtract centre of mass velocities (7-velocities) for each 
source to account for differences between the data sets due 
to the various telescopes, instruments, choice of lines and 
methods of radial velocity determination used. The radial 
velocities of a total of 48 different lines were used in the 
calculation of our ephemeris - the majority of these can be 
found in Table 2 of Gies & Bolton (1982); we also use addi- 
tional He 1 lines (A4921, A5015, A5047, A5875, A6678), oxy- 
gen lines (OnA4349, OnA4366, OmA4650), MgnA4481 
and Nll,A4630. As with previous authors, all lines were 
treated equally, although we note that inhomogeneities in 
the atmosphere of the supergiant (due to the variety of ve- 
locities in the stellar wind) cause greater shifting in the red 
Hcl lines than the blue - unfortunately it is not possible 
to correct for this without knowing exactly which authors 
used which lines for each radial velocity measurement. 

With so many radial velocities from fourteen different 
sources we were then able to omit those spectra for which 
O-C > 20 km/s and still have a total of 421 points from 
which to construct a radial velocity curve. We have also 
omitted the results of Seyfert & Popper 1941 for reasons 
discussed in LaSala et al. (1998). The resultant radial ve- 
locity curve is shown in Figure 1; residuals are also plot- 
ted and show a definite lack of structure. Weighting all 
points by the inverse of their r.m.s. error and assuming 
an orbit of zero eccentricity the orbital parameters were 
obtained using the FOTEL3 software. They can be found 
in Table 3, along with the elements calculated using the 
photometric data. We calculate the orbital period to be 
Porb,s P = 5.599829 ± 0.000016 and define T as the time 
of superior conjunction of the black hole. 

Our new, high quality spectra have enabled us to ex- 
tend the baseline for the ephemeris by over 400 days. As 
well as improving the accuracy of the radial velocity curve 
it also emphasizes the stability of the orbit of Cygnus X-l 
over long time intervals. 



Table 1. Radial velocities for Crimean observations of 
1997 



JDh (2400000+) V r (km/s) r.m.s. error 



50615.4956 


-48.5 


1.0 


50623.4675 


66.4 


1.4 


50624.4585 


-14.3 


1.2 


50625.4553 


-63.3 


2.3 


50626.4792 


-72.7 


1.3 


50648.4623 


-78.5 


1.3 
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50650.5054 


56.8 


1.3 


50651.5033 


62.9 


1.4 


50653.5095 


-65.4 


1.1 


50658.4643 


-26.1 


1.0 


50660.3914 


-52.1 


1.6 


50661.3882 


33.5 


1.7 


50663.4452 


4.2 


1.0 


50667.4747 


77.3 


1.6 


50668.4164 


66.2 


1.7 


50669.4382 


10.0 


2.5 


50675.3550 


-44.7 


1.2 


50676.4727 


-86.0 


1.4 


50677.4751 


-28.8 


1.4 



2.2. Photometry 

The photometric UBV data were obtained on the 60cm 
telescope at the Crimean Laboratory (Nauchny, Crimea) 
of the Sternberg Astronomical Institute during 1971-1997. 
The pulse-counting photometer was used except for eight 
nights in July- August 1997 when it was changed to a CCD 
photometer with ST6V Camera (SBIG). All observations 
were made with reference to the standard BD +35°3816, 
V = 9 m .976,5 - V = m .590,[/ - B = m .064 (Lyuty 
1972). We therefore have a homogeneous photometric data 
set containing more than 800 observtions. For this reason 
we did not use the photometric data of other authors ex- 
cept those of Khaliullin & Khaliullina (1981) which were 
obtained with the same telescope, photometer and local 
standard. 

Using a DFT we obtained the photometric period in 
each of the U, B and V filters. The mean for the three 
filters, P or b,phot — 5.599836 ± 0.000024 days, is in remark- 
able agreement with the spectroscopic period. We do not 
present the power spectrum here as the peak correspond- 
ing to the photometric period (and its aliases) is the only 
significant peak in the frequency range from zero up to 
the Nyquist frequency 

Previously Kemp et al. (1987) gave the most accurate 
value for the photometric period and the epoch of primary 
minimum (superior conjunction of X-ray source): 
Min I = JD2441163.631(±0.005) + 5.59985(±0.00012). 
The low accuracy of the period value may be due to 
combining different data sets. However, using the same 
data Lloyd & Walker (1989) obtained P orb = 5.59982 ± 
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TD ("2400000+) 


Sourc6 


No. of spectra 


^y_ vpIopi tv (km /s ) 


O-C (km /si 


41 159-41588 

t:J. J.tJ^/ j: luUU 


Smith pf al 1Q7^ 


g 


0.1±2.2 


6.3 


41 1 63-41 955 


Webster Mnrdin 1 Q79 


16 


0.1±2.7 


10.5 


41213-44795 


Gies, Bolton, 1982 


78 


-1.8±0.8 


6.6 


41214-41477 


Brucato, Kristian, 1973 


12 


1.1±1.7 


5.8 


41269-41012 


Mason et al., 1974 


14 


2.5±2.3 


8.7 


41515-42670 


Walker, Yang, Glaspey, 1978 


13 


-0.6±1.0 


3.7 


41844-41290 


Brucato, Zappala, 1974 


17 


-3.9±1.8 


7.6 


42205-42910 


Abt, Hitzen, Levi, 1977 


79 


-7.9±1.1 


7.9 


43090-44768 


Aab, 1983b 


24 


-6.3±1.3 


5.9 


44513-45895 


Ninkov et al., 1987 


84 


-5.6±0.7 


6.6 


46332-46635 


Sowers et al., 1998 


14 


-9.8±1.4 


5.1 


49217-49538 


Canalizo et al., 1995 


6 


10.6±1.9 


3.6 


50228-50255 


LaSala et al., 1998 


33 


-0.9±1.8 


10.2 


50615-50677 


Crimea 1997 (this work) 


20 


-4.1±1.9 


8.4 



0.00005, which coincides with the more accurate value 
of 5. 59982 (±.00004) (Voloshina, Lyuty & Tarasov 1997). 
Nadzhip et al. (1996) also improved the epoch of primary 
minimum to 
JD 2441163. 547±0.005. 

Our new epoch of superior conjunction (Table 3) cor- 
responds to primary minimum at JD2441 163.529 ± 0.009, 
so, we can give the most accurate ephemeris as: 
Min I = JD2441163.529(±0.009)+5.599829(±0.000016)£; 

Using this ephemeris we have constructed the mean 
UBV curves for our 27 years' worth of photometric data 
(Figure 2). To increase the homogeneity of the data set 
we have not used all available measurements, but only 1-5 
per night, totalling 827 UBV measurements. The observed 
mean curves are compared with the theoretical ellipsoidal 
curves of Balog, Goncharskij & Cherepashchuk (1981), for 
which i = 50°, q = 0.33, fj, = 0.9, T* = 30000if. 

The discrepancy between the observed mean light 
curves and the theoretical curve is probably due to the 
stellar wind. Were the accretion disk to be responsible 
then we would expect to see the maximum additional 
emission in the U-band (assuming disk temperature of 
~ 10000 K) at inferior conjunction of the black hole (phase 
0.5). Clearly this is not the case, the maximum additional 
emission being in the [/-band at phase 0.05 and 0.5. Fur- 
ther analysis of the photometry is in preparation. 

3. Conclusion 

We have refined the orbital period of Cygnus X-l and 
remain within agreement of the results of many other au- Fi §- 2 - Observed mean lightcurves for 1971-1997 com- 
thors. Our value of 5.599829 days for the orbital period P ared with theoretical ellipsoidal curve 
also provides an accurate fit to 27 years' worth of pho- 
tometry. Further mutiwavelength analysis will follow in 
another paper. 
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Fig. 1. Top: Radial velocity curve using accumulated data compared with theoretical curve from Balog, Goncharskij 
& Cherepashchuk (1981), Bottom: Residuals 



Table 3. The Orbital Elements of Cygnus X-l 



Element 


Spectroscopic 


Photometric 


Period (days) 


5.599829±0.000016 


5.599836±0.000037 


T (JDh) 


2441874.707±0.009 


2441163.529±0.009 


K x (km/s) 


74.93±0.56 




/(M)(M ) 


0.244 





ity of this work was accomplished on the Sussex STARLINK 
node. 
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